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We  study  the  following  problem:  given  a  strongly  connected  digraph,  ^d 
minimal  strongly  connected  spanning  subgraph  of  iL  Our  main  result  is  a  ppalle! 
algorithm  for  this  problem,  which  runs  in  polylog  parallel  time  and  uses  O  (n*)  pro 
cessors  on  a  PRAM.  Our  algorithm  is  simple  and  the  major  tool  it  uses  is  comput 
ing  a  minimum-weight  branching  with  zero-one  weights.  We  also  present  sequen 
tial  algorithms  for  the  problem  that  run  in  time  O  (m  +n  logn). 
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1.  Introduction 

The  transitive  reduction  problem  for  soongly  connected  digraphs  is:  given  a  strongly  con¬ 
nected  digraph  G,  hnd  a  minimal  strongly  connected  spanning  subgraph  of  it,  i.e.,  a  strongly  con¬ 
nected  spanning  subgraph  for  which  the  removal  of  any  arc  destroys  strong  connectivity.  We  are 
looking  for  a  minimal  subgraph  because  the  problem  of  finding  a  minimum  subgraph  with  the 
same  transitive  closure  is  A/P-hard  [GJ]. 

There  is  an  obvious  sequential  algorithm  for  solving  this  problem:  scan  the  arcs  one  by 
one;  at  step  i  test  if  the  i-th  arc  can  be  removed  without  destroying  strong  connectivity.  If  so, 
remove  it  and  update  the  digraph.  This  algorithm  has  complexity  0((n+m)^),  where  n  is  the 
number  of  vertices  of  the  input  graph  and  m  is  the  number  of  arcs.  A  simple  modification  is  to 
initially  reduce  the  number  of  arcs  to  at  most  2n-2  by  taking  the  union  of  a  forward  and  an 
inverse  branching  (defined  below).  This  reduces  the  running  time  to  0  (n^). 

The  problem  studied  here  is  reminiscent  of  the  well-studied  problem  of  finding  a  maximal 
independent  set  of  vertices  in  a  graph,  for  which  several  parallel  algorithms  have  appeared  in  the 
literature  ([KW],[Lu],[ABI],[GS]).  Two  common  features  are  that  there  is  a  simple  sequential 
algorithm  for  it  that  seems  hard  to  parallelize  and  that  the  related  optimization  problem 
(minimum  vs.  minimal)  is  NP-hard. 

We  can  define  the  following  independence  relation  on  the  arcs  of  a  strongly  connected 
digraph,  G:  a  set  of  arcs  is  independent  if  it  can  be  removed  without  destroying  strong  connec¬ 
tivity  of  G.  Using  this  definition,  finding  a  transitive  reduction  of  G  is  equivalent  to  removing  a 
maximal  independent  sets  of  arcs  from  G.  A  property  that  sets  our  problem  apan  from  the  maxi¬ 
mal  independent  set  problem  is  that  in  our  case  independence  of  a  set  is  not  guaranteed  when 
every  pair  of  elements  in  it  is  independent 

Our  problem  can  be  expressed  as  the  determination  of  a  maximal  independent  set  in  an 
independence  system  as  defined  by  Karp,  Upfal  and  Wigderson  ([KUW]).  The  problem  com¬ 
puted  by  a  "rank  oracle"  in  this  case  is  ^P-hard,  but  an  "independence  oracle"  is  easy  to  compute 
in  NC.  Following  the  method  described  in  [KUW]  this  automatically  yields  a  randonmed  paral¬ 
lel  algorithm  that  uses  a  polynomial  number  of  processors  and  runs  in  time  0(Vn  -log^n)  (for 
some  constam  c). 

In  this  paper  we  present  parallel  and  sequential  algorithms  for  this  problem.  Our  first  par  al¬ 
lel  algorithm  runs  in  time  OClog^n)  and  uses  O(n^)  processors  on  a  GIEW  PRAM.  We  then 
presem  an  improved  implementation  of  one  of  the  steps  in  the  algorithm  that  leads  to  a  parallel 
algorithm  that  runs  in  0(los*n)  time  v/ith  the  same  processor  bound.  Both  of  these  algorithms 
can  be  speeded  up  bv  a  logn  factor  if  we  use  a  QICW  PRaM;  we  assume  here  the  COMMON 
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concurrent-whte  model  in  which  all  processors  participating  in  a  concurrent  write  must  write  the 
same  value  [KR],  The  processor  bound  of  O(n^)  represents  the  number  of  processors  needed  to 
multiply  two  n  by  n  matrices  in  0  Oogn)  time  on  a  CREW  PRAM  by  the  straightforward  parallel 
matrix  multiplication  algorithm.  It  is  possible  that  the  processor  bound  can  be  improved  by  using 
sophisticated  techniques  for  multiplying  n  by  n  matrices  (see  e.g.,  [CW]);  we  do  not  elaborate  on 
this. 

The  major  tool  that  our  algorithms  use  is  computing  a  minimum-weight  branching  with 
zero-one  weights.  Central  to  our  algorithms  is  a  proof  that  two  suitable  applications  of  this  tool 
are  guaranteed  to  reduce  by  half  the  number  of  arcs  still  to  be  removed.  We  also  present  two 
sequential  algorithms  for  the  problem,  each  of  which  runs  in  time  Oim+n  \ogn).  This  is  an 
improvement  over  the  straightforward  algorithm  mentioned  above. 

The  transitive  reduction  problem  is,  in  some  sense,  a  dual  of  the  minimum  strong  augmen¬ 
tation  problem  -  add  a  minimum  set  of  arcs  to  a  digraph  to  make  it  strongly  connected.  A  linear 
time  sequendal  algorithm  was  given  for  this  problem  by  Eswaran  and  Taijan  ([ET]),  and  a  paral¬ 
lel  algorithm  running  in  OQogn)  time  with  0(n^)  processors  on  a  CRCW  PRAM  was  given  by 
Soroker  ([So]). 

Our  problem  extends  naturally  to  general  digraphs:  given  a  digraph  G,  find  a  minimal  span¬ 
ning  subgraph  of  it  whose  transitive  closure  is  the  same  as  that  of  G.  A  sequential  algorithm  for 
this  problem  in  the  case  that  G  is  acyclic  is  given  in  [AGU]  and  can  be  parallelized  in  a  straight¬ 
forward  manner.  Combining  it  with  our  algorithms  we  obtain  parallel  algorithms  (with  the  same 
complexities  as  stated  above)  for  the  transitive  reduction  problem  on  general  digraphs.  We  point 
out  that  these  parallel  algorithms  are  good  with  respect  to  the  state  of  the  art,  since  the  problem 
solved  is  at  least  as  hard  as  testing  reachability  from  one  vertex  to  another  in  a  digraph,  and  the 
best  NC  algorithm  currently  known  for  this  requires  on  the  order  of  Af(n)  processors,  where 
M(n)  is  the  number  of  processors  needed  to  multiply  two  nby  n  Boolean  matrices  in  0(\ogn) 
time. 

We  note  that  the  name  "transitive  reduction"  was  given  to  a  different  problem  by  Aho, 
Garey  and  Ullman  ([AGU]),  and  should  not  be  confused  with  our  definitiort  Given  a  digraph  G, 
they  ask  for  a  digraph  with  a  minimum  number  of  arcs  (not  necessarily  a  subgraph  of  G)  whose 
U«uu>iuvc  c!?‘:’jre  is  the  same  as  that  of  G.  This  definition  agrees  with  curs  when  G  is  acyclic 

Definitions 

Let  G  be  a  strongly  cormected  digraph.  A  forward  (inverse)  branching  rooted  at  x  is  a 
spanning  tree  of  G  in  which  x  has  in-degree  (out-degree)  zero  and  all  other  vertices  have  in¬ 
degree  (out-degree)  one.  A  branching  is  either  a  forward  or  an  inverse  branching.  Throughout 
this  paper  the  root,  x,  will  be  some  (arbitrarily)  fixed  vertex  of  the  input  digraph,  and  the  set  of  aU 
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branchings  will  be  taken  to  be  only  those  rooted  at  x. 

An  arc,  e,  is  G -redundant  (or  simply  redundant  when  the  graph  is  clear)  if  G-le)  is 
strongly  connected.  Arc  e  is  G -essential  (or  essential )  if  it  is  not  redundant.  Let  //  be  a  sub¬ 
graph  of  G.  Let  rn(H)  denote  the  number  of  G- redundant  arcs  in  H.  When  H=Gwc  will  use  the 
shorthand  r(G). 

An  H-philic  ( H -phobic )  branching  in  G  is  one  that  has  the  greatest  (smallest)  number  of 
arcs  in  common  with  H  over  all  branchings  (rooted  at  x)  in  G. 

Our  model  of  parallel  computation  is  the  Parallel  Random  Access  Machine  (PRAM),  which 
consists  of  a  collection  of  independent  processing  elements  commumcating  through  a  shared 
memory.  For  a  survey  on  the  PRAM  model  and  PRAM  algorithms  see  [KR]. 


2.  The  Transitive  Reduction  Algorithm 

Our  basic  algorithm  is  based  solely  on  computing  philic  and  phobic  branchings.  The  fol¬ 
lowing  lemma  explains  how  these  branchings  are  computed: 

Lemma  0;  An  //-philic  (//-phobic)  branching  can  be  computed  by  a  minimum -weight  branching 
computation  with  zero-one  weights. 

Proof:  Assign  weight  0  (1)  to  every  arc  in  //  and  weight  1  (0)  to  all  other  arcs.  [] 

Such  a  minimum-weight  branching  can  be  computed  in  time  0(\o^n)  using  G(n^)  proces¬ 
sors  on  a  ORCW  PRAM  by  Lovasz’s  method  ([Lo]).  On  a  CREW  PRAM,  this  algorithm  runs  in 
O  Oog^n)  time. 

Proposition  1:  An  arc  of  G  is  essential  if  and  only  if  it  is  the  unique  arc  crossing  some  directed 
cut  of  G. 

Proposition  2:  The  union  of  a  forward  branching  and  an  inverse  branching  of  G  is  a  strongly 
connected  spanning  subgraph  of  G. 

Proposition  3:  Let  G'  be  a  strongly  connected  spanning  subgraph  of  G.  Then  e  is  G'-redundant 
only  if  it  is  G-redundam. 

Lemma  1:  Let  F  be  a  forward  branching  in  G  and  let  /  be  an  F-philic  inverse  branching  in  G. 
Let  G'=Ft^/ 'Ihcn  the  arcs  of /-F  ate  all  G'-essential. 
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Proof:  Let  eel-F.  Assume  G'-[e)  contains  some  inverse  branching,  /'.  Then  /'  has  one  more 
arc  in  common  with  F  than  /  does  (since  all  branchings  have  the  same  number  of  arcs).  But  this 
contradicts  the  fact  that  /  is  F-philic.  Thus  G'-[e}  contains  no  inverse  branching  and  is  therefore 
not  strongly  conneaed.  [] 

A  cut  leaving  S  is  the  set  of  arcs  extending  from  S  to  Y.(G)-~S  in  a  digraph,  G,  and  its  cardi¬ 
nality  is  denoted  by  5n(S). 

Theorem  1  (Edmonds’  Branching  Theorem  ([Ed])): 

Ut 

k  =min  {5c(S)  1  xeS  ,  S^V(G) }. 

Then  G  contains  k  arc-disjoint  forward  branchings  (rooted  at  x). 

Lemma  2:  For  every  strongly  connected  digraph,  G,  there  exists  a  forward  branching,  F,  of  G 
such  that  rQ(F)^'Ar  (G). 

Proof:  Let  G'  be  obtained  from  G  by  duplicating  all  essential  arcs.  Let  S  be  a  proper  subset  of 
V(G)  containing  x.  We  claim  that  5(;'(S)>2.  This  is  because  the  cut  leaving  S  must  contain  at 
least  one  duplicated  essential  arc  of  G  or  at  least  two  redundant  arcs  (by  proposition  1).  There¬ 
fore,  by  theorem  1,  there  are  two  arc-disjoint  forward  branchings  in  G'  (each  corresponding  to  a 
branching  in  G),  one  of  which  must  contain  at  most  half  of  the  (unduplicated)  G-redundant  arcs.[] 

Theorem  2:  Let  F  be  the  set  of  redundant  arcs  in  G.  Let  F  be  an  F-phobic  forward  branching 
and  let  /  be  an  F-philic  inverse  branching.  Let  G'=F\jl.  Then  r  (G')^'/^(C/). 

Proof:  First  note  that  by  proposition  2,  G'  is  strongly  connected.  By  lemma  2  and  proposition  3, 
rG’(F)^GiF)^^^(P)-  By  lemma  l.r(G>/-G-(F).  Therefore  r(G')S'/4r(G).  [] 

It  is  an  immediate  consequence  of  theorem  2  that  the  following  NC  algorithm  gives  a  transitive 
reduction  of  G: 

Repeat 

(1)  F  <-  set  of  redundant  arcs  in  G 

(2)  F  <-  F-phobic  forward  branching  in  G 

(3)  I  <-  F-philic  inverse  branching  in  G 

(4) G4-Fu/ 
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until  R=^ 

(5)  output  G  (it  is  a  transitive  reduction  of  the  input  digraph) 

By  Theorem  2  the  repeat  loop  runs  O  Gogn)  times,  where  n  is  the  number  of  vertices  in  G.  Steps 
(2)  and  (3)  are  implemented  with  Lovasz’s  minimum-weight  branching  algorithm  Gemma  0). 
The  straightforward  implementation  of  step  (1)  is  to  perform  a  strong  connectivity  test  (transitive 
closure)  with  each  vertex  of  the  graph  deleted  in  turn,  which  requires  n-Min)  processors.  In  the 
next  section  we  shall  show  how  to  perform  this  step  more  efficiently. 


3.  Efficient  Classification  of  Arcs 

In  this  section  we  give  parallel  algorithms  to  classify  the  arcs  of  G  as  essential  or  redundant 
in  poly-log  time  using  only  0(n^)  processors.  In  section  3.1  we  provide  a  simple  polylog  time 
parallel  algorithm  using  O(n^)  processors.  In  section  3.2  we  provide  a  faster  algorithm  using  tree 
contraction  [MR]. 

3.1.  Finding  Redundant  Arcs  Using  Minimum  Weight  Branchings 

Let  ^  ( )  be  the  set  of  essential  arcs  contained  in  all  forward  (inverse)  branchings.  It  fol¬ 
lows  from  proposition  2  that: 

Proposition  4:  An  arc  is  essential  if  it  is  either  in  Ej  or  in  Ei  (or  both). 

Lemma  3:  Let  //  be  a  set  of  arcs  containing  Ef  and  let  F  be  an  //-phobic  forward  branching  in  G. 
Then  l(Fp^//)-£/l  ^  \H-Ef\. 

Proof:  Let  G'  be  obtained  from  G  by  duplicating  all  the  arcs  in  Ef.  As  in  lemma  2,  there  exist 
two  arc-disjoint  forward  branchings  in  G'  (corresponding  to  branchings  in  G).  one  of  which  con¬ 
tains  at  most  half  the  arcs  of //-£y.  [] 

Therefore  Ef  (and  similarly  £,)  can  be  computed  by  the  following  algorithm: 

(l)//<-G 

repeat  steps  (2)  and  (3)  \  Ig  nii  times 
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(2)  F  <-  //-phobic  forward  branching  in  G 
(4)  output  H  (this  is  the  set  E/) 


This  algorithm  requires  logn  applications  of  Lovasz’s  minimum  weight  branching  algo¬ 
rithm.  which  runs  in  0  Gog^/i)  parallel  time  on  a  CRCW  PRAM  with  0  (n^)  processors.  Thus  we 
can  use  this  algorithm  to  find  all  redundant  arcs  in  OOog^n)  parallel  time  on  a  OiCTW  PRAM 
with  O(n^)  processors.  This  in  turn  leads  to  a  transitive  reduction  algorithm  that  runs  in 
O  Gog^n)  parallel  time  on  a  CRCW  PRAM  with  O  (n^)  processors. 


32.  Finding  Redundant  Arcs  Using  Tree  Contraction 

Let  r  be  a  fixed  root  of  a  directed  graph  G=(V',£).  We  call  arc  (v.w)  an  out-bridee  if  (v.w) 
is  on  every  path  from  r  to  w,  and  an  in-bridee  if  (v.w)  is  on  every  path  from  v  to  r.  Let  G  be  the 
set  of  out-bridges  of  G,  and  I  the  set  of  in-bridges  of  G.  Then  the  set  of  redundant  arcs  is  the  set 
E-UkjO). 

Let  R  be  a  forward  branching  rooted  at  r.  Then  every  out-bridge  of  G  lies  in  B.  We  can  view 
5  as  a  rooted  directed  tree  B=(y,B',r).  For  a  vertex  v  in  V-  (r),  we  denote  by  parent iy),  the 
parent  of  v  in  fl.  A  vertex  v  is  active  if  there  is  a  path  from  r  to  v  that  avoids  arc  (parent (y),v). 
Similarly,  a  non-tr';e  arc  (w.v)  is  active  if  it  lies  on  a  path  from  r  to  v  that  avoids  arc 
(parent  (v),v). 

Lemma  4:  Let  B=(y,E',r)  be  a  forward  branching  in  a  directed  graph  G.  A  tree  arc 
e=(parent(v),v)  in  B  is  an  out-bridge  of  G  if  and  only  if  v  is  not  active. 

Proof:  If  e  is  an  out-bridge  of  G  then  every  path  from  r  to  v  passes  through  e.  Thus  v  cannot  be 
active.  Conversely,  if  e  is  not  an  out-bridge,  then  there  exists  a  path  from  r  to  v  that  avoids  e  and 
hence  v  must  be  active.[] 

We  now  give  an  algorithm  to  identify  all  active  vertices,  and  hence  all  out-bridges,  using 
tree  contraction  [MR].  An  analogous  computation  on  an  inverse  branching  rooted  u  r  gives  the 
in-bridges,  from  which  we  can  compute  the  redundant  arcs  in  G. 

We  shall  use  a  variant  of  uce  contraction  proposed  in  in  which  the  basic  operation  is 
shrink,  which  we  now  define.  A  leaf  chain  in  a  rooted  tree  T Hy,E,r)  is  a  path  <v  j ,  •  •  • ,  v/>  such 
that  each  v,,  i  >  1  has  exactly  one  incomirig  arc  and  one  outgoing  arc  in  T.  v  i  has  either  no  incom¬ 
ing  arc  or  more  than  one  outgoing  arc  in  T,  and  v/  is  a  leaf  in  T.  We  will  call  v  |  the  root,  and  v/ 
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the  leaf  of  the  leaf  chain.  Note  that  every  leaf  in  7  is  part  of  a  leaf  chain,  possibly  a  degenerate 
one  (if  1=2). 

The  shrink  operation  applied  to  a  rooted  tree  T  =(V  E,r)  removes  all  vertices  in  each  leaf 
chain  in  T  except  the  root  of  the  leaf  chain.  It  can  be  shown  that  O  Gog«)  applications  of  the 
shrink  operation  suffice  to  reduce  any  /i-node  tree  to  a  single  node  [Ra]. 

We  now  develop  an  algoritlim  ShrinkfP)  for  identifying  out-bridges  for  the  case  when  the 
forward  branching  is  a  simple  path.  We  shall  then  use  this  to  find  the  out-bridges  in  leaf  chains 
while  implementing  the  shrink  operation  in  a  tree  contraction  algorithm  to  find  out-bridges  in  G 
given  an  arbitrary  forward  branching. 

The  input  to  algorithm  Shrink(P)  is  a  directed  graph  ?=(V,£)  consisting  of  a  directed  path 
p=<l,2,  •••,:>,  together  with  a  collection  of  forward  arcs  of  the  form  <j.  and  a  collec¬ 
tion  of  back  arcs  of  the  form  The  algorithm  Shrink(P)  will  identify  all  active  vertices, 

thereby  giving  the  out-bridges  in  p.  Note  that  P  is  allowed  to  have  two  arcs  of  the  form 
one  of  which  is  a  forward  arc  and  the  other  lies  in  p.  We  will  need  this  when  we  apply  algorithm 
Shrink(P)  to  the  general  problem  of  finding  out-bridges  in  a  graph  with  an  arbitrary  forward 
braiKhing. 

We  now  make  a  series  of  observations. 

Observation  1:  Every  forward  arc  is  active. 

Let  p(u)  be  the  subgraf^  of  G  induced  by  vertices  u  through  t.  For  each  vertex  v  in  p(a), 
let  V— if  u  is  reachable  from  v  inp  («).  Let  reach  («)  be  the  set  of  vertices  v  in  p  (u)  with  v-*u. 

Observation  2:  Reach  (u)  is  a  single  interval  of  the  form  [u,u'].  Farther  a  vertex  v*u  is  in 
reachiu)  if  and  only  if  there  exists  a  sequence  of  back  arcs  &,=(u,,v, ),/=!,  •  •  ,k  such  that  vi=w, 

and  tt, ^,+1,1=1,  •  •  •  ,k-\. 

Lemma  5:  A  vertex  u  is  active  if  and  only  if  there  is  a  forward  arc  (.k,l)  with  k<u  an6  I  in 
reach  (w). 

Proof:  Let  u  be  an  active  vertex.  Then  there  is  a  path  q  from  the  root  to  u  that  avoids  arc  (u-1  ,u). 
This  in  turn  implies  that  q  must  contain  a  forward  arc /  =(k,l)  with  k  <u,  Itu  and  with  u  reachable 
from  /  using  only  arcs  in p ,  i).  Hence  /  must  be  in  reach (u). 

Conversely  suppose  there  is  a  forward  arc  /  =ik,l)  with  k  <u  and  I  in  reach  (u).  Hence  there 
is  a  path  q  from  /  to  u  using  only  arcs  in  p  (u).  Then  the  path  consisting  of  arcs  in  p  from  the  root 
to  k,  followed  by  arc / and  then  the  path  qisi  path  from  1  to  u  that  avoids  arc  (u-l,u).  Hence  u 
must  be  an  active  vertex.  [] 

Observations  1  and  2  and  Lemma  5  together  give  us  the  following  algorithm  to  find  all  out- 
bridges  when  the  forward  branching  is  a  simple  path. 


Shrink(P); 

1.  Find  reach  («)  for  each  vertex  u  as  follows: 

a)  For  each  back  arc  b=(j,;)  find  a  back  arc  next{b)=(j^ f)  with  /  in  [J,i  ]  and  maximum  If 
i'^  then  set  next(b)=^. 

b)  Form  an  auxiliary  graph  with  a  vertex  for  each  back  arc  b  and  an  arc  from  b  to  next  (b),  if 
nextib)  exists.  This  auxiliary  graph  is  a  forest  of  trees. 

c)  For  each  vertex  u,  pick  some  back  arc  b=(y,u)  incident  on  u,  and  find  the  root  b'  of  the 
tree  it  belongs  to.  Let  6'  be  the  back  arc  (x.y).  Set  reach  («>=[tt,x  ]. 

If  there  is  no  back  arc  incident  on  u  set  reach  (u)=^u,u]. 

2.  For  each  vertex  u,  find  a  forward  arc  /  with  I  in  reachiu)  and  with  minimum  k.  If  k<u 
mark  u  as  active. 

3.  For  each  vertex  u  that  is  not  active,  mark  (p  (u),u)  as  an  out-bridge. 

We  now  show  how  to  implement  each  of  the  steps  in  the  algorithm  efficiently  in  parallel. 
Step  3  can  be  implemented  trivially  in  constant  time  with  t  processors.  The  following  method 
implements  step  2  in  O  Gogt)  time  with  a  number  of  processors  linear  in  the  size  of  P\  Initially 
we  determine,  for  each  vertex  u,  the  forward  arc  <v,u)  with  minimum  v  (if  such  an  arc  exists).  It 
is  straightforward  to  compute  tliis  in  O  Gog/)  time  with  a  linear  number  of  processors.  Then  by  a 
doubling  computation  we  compute,  for  each  interval  Iu,m+2''],  ^og^  ,  the  forward 

arc  (v,x)  with  minimum  v  such  that  x  is  in  the  interval  [u,u-f2-'].  This  computation  can  oe  done  in 
O  Gog/)  time  with  a  linear  number  of  processors  on  a  CREW  PRAM.  Any  interval  [i,j  ],  1^'  <j^ 
can  be  written  as  the  overlapping  union  of  two  of  the  previously  computed  intervals,  and  hence 
each  vertex  can  now  find  a  forward  arc  as  required  in  step  2  in  constant  time. 

Step  1  can  be  implemented  in  O  Gog/)  time  with  a  liner'  number  of  processors  on  a  CREW 
PRAM  as  follows.  Step  la  can  be  performed  in  a  manner  analogous  to  step  2.  Step  lb  can  be 
implemented  in  constant  time  with  a  linear  number  of  processors.  Step  Ic  can  be  implemented  by 
pointer  jumping  in  O  Gog/)  time  with  a  linear  number  of  processors.  Thus  we  have  a  parallel 
algorithm  for  Shiink(P)  that  runs  in  O  Gog/)  time  with  a  linear  number  of  processors  on  a  CREW 
PRAM. 

We  now  incorporate  the  Shrink  algorithm  in  the  following  tree  contraction  algorithm  that 
finds  the  out-bridges  in  an  arbitrary  forward  branching  of  a  directed  graph  G  rooted  at  r.  The 
algorithm  constructs  a  sequence  of  pairs  (Gt.Tit),  where  Gt  is  a  digraph  and  T^t  is  a  forward 
branching:  G\  is  the  input  digraph  and  Ti  is  a  forward  branching  of  G  rooted  at  a  fixed  vertex  r. 
Iteration  k  identifies  the  leaf  chains  of  determines  the  out-bridges  of  G^  within  those  leaf 
chains,  deletes  all  the  vertices  of  the  leaf  chains  except  their  roots,  and  then  performs  a  transitive 
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closure  computation  and  adds  appropriate  arcs  to  ensure  that  the  out-bridges  in  are  precisely 
the  out-bridges  of  G  not  yet  identified. 

Outbridges(G=(V£,r),T); 

Input:  A  directed  graph  G  ={V,E)  with  a  forward  branching  T rooted  at  r,  IV  l=n. 

Repeat 

1.  Find  out-bridges  in  the  leaf  chains  ofT: 

For  each  leaf  chain  I  in  T  pardo 

Lei  r  be  the  root  of  /  and  r'  the  leaf  of  1.  Let  L'  be  the  subgraph  of  G  induced  by  vertices  in  /. 

a)  Form  L  from  L'  by  introducing  a  forward  arc  (r,y)  for  each  non-tree  arc  (x,y)  in  G  with  y 
in  V(l)-  {/}  andxnot  in  V(l). 

b)  Apply  ShrinJdL)  to  find  the  out-bridges  in  L  and  label  these  as  out-bridges  of  G. 

2.  Remove  leaf  chains  from  T: 

a)  Form  tht  graph  H  with  vertex  set  V  and  arc  set  the  arcs  in  all  leaf  chains  and  all  non-tree 
arcs  of  G. 

b)  Form  M,  the  adjacency  matrix  of  H,  and  form  the  transitive  closure  Af  *  of  Af. 

c)  For  each  vertex  v,  determine,  using  A/*,  the  set  of  vertices  from  which  v  is  reachable  in 
H.  For  each  such  vertex  w,  introduce  an  arc  (h'.v)  in  G. 

d)  For  each  vertex  t  that  is  the  head  of  some  leaf  chtin,  delete  all  incoming  non-tree  arcs  to 
proper  descendants  of  t.  Collapse  ail  of  these  proper  descendants  into  t.  Delete  any  self¬ 
loops  in  this  graph. 

until  T  =4 


Generalizing  our  earlier  notation  for  the  case  when  the  forward  branching  is  a  simnle  oath, 
we  now  let  p(u)  be  the  subgraph  of  G  induced  by  those  vertices  that  lie  in  the  subtree  of  T  rooted 
at  u.  For  each  vertex  v  in  p(u),  let  v-»u  if  u  is  reachable  from  v  inp(u).  Let  reach  (u)  be  the  set 
of  vertices  v  in  p(u)  with  v-*u. 

The  following  lemma  is  a  straightforward  generalization  of  Observation  2  and  Lemma  5 
(here  a  vertex  v  is  a  descendant  of  a  vertex  u  if  u=v  or  if  there  is  a  directed  path  from  u  to  v  in  T; 
otherwise  v  is  a  non-descendant  of  u). 

Lemma  6:  A  vertex  u  in  G  is  active  if  and  only  if  there  is  an  arc  (x.y)  with  x  a  non-descendant  of 
u  and  with  y  in  reach  (u). 
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Let  G  be  a  dire  xd  graph  with  a  forward  branching  T  rooted  at  r,  and  let  v  be  a  vertex  in  G. 
An  active  oath  to  v  is  a  path  p  from  r".  v  consisting  of  an  iriiial  path  p'  using  tree  arcs  from  r  to 
a  non- descendant  z  of  v  followed  by  an  intermediate  path  consisting  of  a  single  non-tree  arc  a 
from  z  to  a  descendant  y  of  v  followed  by  a  final  path  p"  from  y  to  v  using  only  arcs  connecting 
descendants  of  v. 

Observation  3:  Vertex  v  is  active  if  and  only  if  there  is  an  active  path  to  v. 

We  now  prove  some  lemmas  that  will  allow  us  to  establish  the  correemess  of  algorithm 
Outbridges.  As  before  let  G,  and  r,  be  the  graph  and  forward  branching  present  at  the  stan  of  the 
ith  iteration  of  the  repeat  loop  in  the  algorithm;  hence  Gj  and  Ti  are  the  input  graph  together 
with  its  forward  branching,  and  G*  and  T*  are  the  current  graph  and  forward  branching  at  Lhe  start 
of  the  Jkth  iteratioa  Similarly  let  //,  be  the  graph  H  of  step  2a  of  algorithm  Outbridges  constructed 
in  the  ith  iteration  of  the  repeat  loop. 

We  first  note  that  Observation  2  remains  valid  in  each  G*  when  «  is  a  vertex  in  a  leaf  chain 
of  r*.  Wfc  stale  this  in  the  following  observation. 

Observation  4:  Let  u  be  a  vertex  in  a  leaf  chain  /  of  forward  branching  T,  where  for  convenience 
we  assume  that  the  vertices  in  the  leaf  chain  are  numbered  from  1  to  s,  with  1  the  root  of  the  leaf 
chain  and  s  the  leaf  of  the  leaf  chain.  Then  reachiu)  is  a  single  interval  of  the  form  [«,«']. 
Further,  a  vertex  is  in  reachiu)  if  and  only  if  there  exists  a  sequence  of  back  arcs 
••  •  ,k  in  L  (where  L  is  the  subgraph  of  G  induced  by  vertices  in  /)  such  that  vi=u, 
Uk>v,  and  M,^Vi+i,t=l,  •  •  •  ,ik-l. 

Lemma  7:  For  each  k^l,  algorithm  Outbridges  correctly  finds  the  out-bridges  in  the  leaf  chains 
ofG*. 

Proof:  By  Observation  4,  for  a  vertex  u  in  a  leaf  chain  /  of  T*,  reach  («)  in  G*  is  the  same  as 
reach  (u)  in  the  subgraph  of  G^  induced  by  1.  Hence  the  reach  value  of  each  vertex  in  the  leaf 
chain  is  correctly  computed  in  the  Shrink  computation  of  step  lb  in  algorithm  Outbridges. 

By  Lemma  6,  a  vertex  u  in  a  leaf  chain  is  active  if  and  only  if  there  is  an  arc  e=<z,y}  in 
with  z  a  non-descciidant  of  u  and  with  y  in  reach  (u).  Such  an  arc  e  is  either  a  forward  arc  in  the 
leaf  chain  or  is  an  arc  with  z  not  in  the  leaf  chain  and  y  in  the  leaf  chain.  The  former  case  is  the 
same  as  that  used  in  the  Shrink  algorithm.  In  the  latter  case,  (z,y)  will  cause  any  vertex  u  in  the 
leaf  chain  with  y  in  reach  (u)  to  be  active.  Hence  for  the  purpose  of  the  Shrink  algorithm  this  is 
equivalent  to  having  an  arc  from  the  root,  r,  of  the  leaf  chain  to  y.  Thus  the  computation  in  steps 
la  and  lb  of  algorithm  Outbridges  correctly  finds  the  outbridges  in  the  leaf  chains  of  G^.d 

Lemma  8  Let  e=(u,v)  be  an  out-bridge  in  G;i,k>l.  Then  e  is  an  out-bridge  in  G^-i . 

Proof:  First  note  that  if  e  is  an  out-bridge  in  G*,  then  e  lies  in  7*.  Hence  e  lies  in  r*_i,  since 
every  tree  arc  in  7*  is  present  as  a  tree  arc  in  74_i. 
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Suppose  e  is  noi  an  oui-bridge  in  G*_i .  Hence  v  is  an  active  vertex  in  G*_i .  Let  p  be  an 
active  path  to  v  in  G*_i ,  and  let  p  consist  of  an  initial  tree  path  p'  to  a  vertex  x  that  is  a  non¬ 
descendant  of  V,  followed  by  a  non-tree  arc  a={x,y),  where  y  is  a  descendant  of  v,  followed  by  a 
final  path  p"  from  y  to  v  using  only  arcs  connecting  vertices  that  are  descendants  of  v.  We  now 
establish  that  there  must  be  an  active  path  to  v  in  G*.  contradicting  the  assumption  that  e  is  an 
out-bridge  of  G*.  and  thereby  establishing  the  lemma. 

If  p  contains  no  vertex  in  G*_i-G*  then  p  is  an  active  path  to  v  in  G*  as  welL  If  p  contains 
some  vertices  in  Gt-i-G*  then  consider  the  last  vertex  2  on  p  such  that  2  is  in  Gt_i-G*. 

Case  1:  2  is  a  non-descendant  of  v.  Then  2  must  be  x  and  all  vertices  in  p"  lie  in  G*.  Let  r  be  the 
root  of  the  leaf  chain  of  G*_i  to  which  2  belongs.  Then  by  step  2d  of  algorithm  Outbridges,  2  is 
collapsed  into  t  and  hence  the  path  in  G*  consisting  of  the  tree  path  to  r,  followed  by  non-tree  arc 
(r.y),  followed  by  path  p"  is  an  active  path  to  v  in  G*. 

Case  2:  2  is  a  descendant  of  v.  Let  b=(z,a)  be  the  outgoing  arc  from  2  in  p,  and  let  r'  be  the  root 
of  the  leaf  chain  in  G*_i  to  which  2  belongs.  Hence  /'  is  a  descendant  of  a  and  2  is  a  proper  des¬ 
cendant  of  Let  p'"  be  the  portion  ofp"  from  a  to  v.  The  path  p'"  is  a  path  in  G*  as  well. 

Case  2a:  The  vertex  2  is  reachable  from  some  non-descendant  w  of  v  in  //*.  Then  an  arc  (^,2)  is 
introduced  in  step  2c  of  the  algorithm.  If  w  is  in  Gj;  then  the  path  from  r  to  w  followed  by  arc 
(w.a)  followed  by  path  p'"  is  an  active  path  to  v  in  G*.  If  w  is  in  Gi_i-G*  then  the  analysis  of 
Case  1  gives  an  active  path  to  v  in  G*. 

Case  2b:  The  vertex  2  is  not  reachable  from  any  non-descendant  of  v  in  //*.  Now  consider  p". 
This  is  a  path  of  the  form  •  •  •  .«!.*, .vij.  •  •  •  .vj/,.  •  •  •  ,Uc,\,  •  •  •  .«c.*..Vc.i.  •  •  • 

where  the  uij  are  in  G*_i  -G*  and  the  v,  y  are  in  G*.  and  if  y  is  in  G*  the  initial  sequence  of  u  1  ^’s 
is  empty.  All  of  the  ui  j  and  v,^  are  descendants  of  v.  Each  v,,  j  is  reachable  from  ^  ,r  >  1  in 
Ht.  Hence  by  step  2c  of  algorithm  Outbridges.  there  is  an  arc  from  to  j  ,r  >  1  in  G*.  The 

vertex  vi,i  has  an  incoming  arc  from  x  in  G*.  The  remaining  arcs  in  p"  remain  in  G*.  Hence 
there  is  a  path  from  x  to  v  in  G*  that  contains  only  vertices  that  are  descendants  of  v  in  G*.  Hence 
V  is  an  active  vertex  in  G;k-[] 

Lemma  9:  Let  e=(u,v)  be  a  tree  arc  in  Gt,k  >  1  that  is  not  an  out-bridge  in  G*.  Then  e  is  not  an 
out-bridge  in  G*_i. 

Proof:  Since  e  is  not  an  out-bridge  in  G*  there  is  an  active  path  p  to  v  in  G*.  Consider  any  arc 
f=(x,y)  in  p  that  is  not  present  in  G*-i.  If  /was  introduced  in  step  2c  of  algorithm  Outbridges 
then  there  is  a  path  from  x  to  y  in  G*-!  tnat  avoids  all  tree  arcs  in  G*  and  hence  arc  c.  If  /  was 
introduced  in  step  2d  then  there  is  a  path  from  a  descendant  of  x  to  y  in  G^-i  that  avoids  all  tree 
arcs  in  C*.  Hence  there  is  a  path  from  x  to  y  in  G.t_i  that  avoids  arc  e.  Hence  from  p  we  can 
obtain  an  active  path p'  to  v  in  G^-i .  Thus  e  is  not  an  out-bridge  in  Gt-i .[] 
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Lemma  10:  Algorithm  Outbridges  correctly  finds  the  out-bridges  of  G. 

Proof  :  We  show  that  at  the  start  of  each  iteration  of  the  repeat  loop, 

1)  The  out-bridges  identified  so  far  are  ex^-ctly  the  out-bridges  in  the  portion  of  the  input  graph  G 
that  has  been  collapsed  by  the  algorithm. 

2)  An  arc  e  in  the  current  graph  G  is  an  out-bridge  in  this  graph  if  and  only  if  it  is  an  out-bridge  in 
the  original  input  graph. 

The  proof  is  by  inducdon  on  k,  the  number  of  iterations  of  the  repeat  loop. 

Base:  k=\.  The  claim  is  vacuously  true  since  no  out-bridges  have  been  identified  and  the  input 
graph  is  the  same  as  the  current  graph. 

Induction  step:  Assume  that  the  two  claims  are  true  until  the  start  of  iteration  k-\  and  now  con¬ 
sider  the  start  of  iteration  k.  Qaim  1)  follows  by  the  induction  hypothesis  and  t  emma  7.  Claim 
2)  follows  by  the  induction  hypothesis  and  Lemmas  8  and  9.[] 

Finally  we  note  that  algorithm  Outbridges  runs  in  with  0(/i^)  processors  on  a 

CRCW  PRAM.  To  see  the  processor  and  time  bouitds  let  us  analyze  the  time  complexity  of  each 
iteration  of  the  repeat  loop.  By  the  previous  analysis  for  the  time  complexity  of  algorithm  Shrink, 
step  1  runs  in  O  Oogn)  time  with  O  (n^)  processors  on  a  CREW  PRAM.  Steps  2a,  2c,  and  2d  run 
in  0  Qogn)  time  with  0  (n^)  processors  on  a  CREW  PRAM.  Step  2b  runs  in  0  Oogn)  time  with 
O(n^)  processors  on  a  CRCW  PRAM,  and  is  the  most  expensive  step  in  the  repeat  loop.  Since 
the  repeat  loop  is  executed  0  Oogn)  times  we  obtain  the  stated  time  and  processor  bounds  for 
algorithm  Outbridges.  On  a  CREW  PRAM  this  algorithm  runs  in  OGog^n)  time  with  M(n)  pro¬ 
cessors. 

Whether  we  use  a  CREW  model  or  a  CRCW  model  the  time  and  processor  bounds  for 
finding  a  minimum  weight  branching  using  the  algorithm  in  [Lo]  dominate  the  time  and  proces¬ 
sor  bounds  of  algorithm  Outbridges.  Hence  we  can  find  redundant  arcs  within  the  time  and  pro¬ 
cessor  bounds  for  minimum  weight  branchings,  and  thus  the  parallel  transitive  reduction  algo¬ 
rithm  runs  in  0(log^/i)  parallel  time  with  0(/t^)  processors  on  a  CRCW  PRAM  and  in  OGog^n) 
parallel  time  with  the  same  processor  bound  on  a  CREW  PRAM. 


4.  Sequential  Algorithms  for  Transitive  Reduction 

As  in  section  3.2,  let  r  be  a  fixed  root  of  a  directed  graph  G=(V,£),  where  I V  l=n  and 
IE  \  =m.  An  algorithm  for  finding  the  in-  and  out-bridges  is  given  in  [Ta2].  This  algorithm  actu¬ 
ally  does  more:  It  computes  two  forward  branchings  Ti  and  72  having  only  the  out-bridges  in 
common,  and  two  inverse  branchings,  and  T4,  having  only  the  in-bridges  in  common.  This 
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algorithm  can  be  implemented  to  run  in  linear  time  by  using  linear-time  algorithms  for  comput¬ 
ing  nearest  common  ancestors  [HT]  and  maintaining  disjoint  sets  [GT]. 

Let  /?  be  the  set  of  redundant  arcs,  /  the  set  of  in-bridges  and  G  the  set  of  out-bridges.  Hence 
R=E  —  The  following  algorithm  finds  a  transitive  reduction  of  G. 

1.  Pick  a  root  vertex  r  in  G.  Find  a  forward  branching  B  and  an  inverse  branching  B'  in  G  and 
replace  G  by  B  \^B'. 

2.  Repeat 

a)  Construct  two  forward  branchings  T\  and  T2  having  only  the  out-bridges  in  common; 
identify  the  set  of  out-bridges  as  O. 

b)  Construct  two  inverse  branchings  and  T4  having  only  the  in-bridges  in  common; 
identify  the  set  of  in-bridges  as  /. 

c)  Form  the  set  of  redundant  arcs  R^sR  =E-{0\jI). 

d)  For  j=l,2,3,4  form  Si=Ti(^R. 

c)  Choose  Ti  and  Tj  such  that  1:2/ <2.  3<j^  and  has  minimum  cardinality  among 

S\KjS^,S2{jSi,Sn^4,S2(jS4. 

f)  Replace  G  by  Ti\jTj. 
until  R 

The  following  claim  establishes  that  the  repeat  loop  is  executed  only  O  Gogn)  times. 

Lemma  11:  In  step  2e  of  the  algorithm  the  chosen  5,  and  Sj  satisfy  IS,\jSy  I  S(3/4)- 1/?  1 . 

Proof:  For  j=l  ,2,3,4,  let  F,  be  the  set  of  those  arcs  in  Ti  that  arc  not  present  in  any  other  Tj,  and 
let  Fi.3=(-SirV3)VJ^i.  ^2.4=(5'2r>S'4)u^2.^2.3=(52n^3)u^3  and  Fl.4=(5,n‘54)U^4. 
Note  that  P  1.3, P2.4.^2.3  and  P1.4  are  disjoint.  Let  one,  say  P  1.3  be  the  one  of  maximum  cardi¬ 
nality.  Then  we  must  have  lP2,4l+l^2.3  l-*-IPi.4  l^(3/4)  IP  I.  But  .y2L>^4CP2.4tj^2.3t^  j.4. 
which  implies  l52VJ^4l^3/4)’IP  I.  Since  we  also  have  5it^3CPi,3t^.y*i,4t,jP2.3> 
5itj54CPi.3U^i.4Vy*2.4  and  S2t^3C^2.3U^2,4U^i.3.  wc  have  ISn^Ss  1^3/4)- IP  I  if 
P2,4  is  of  maximum  cardinality,  I5it^4l^3/4)’IP  I  if  P2.3  is  of  maximum  cardinality  and 
I  52L>^3  I  ^3/4)- 1 P  I  if  P  1.4  is  of  maximum  cardinality.  Hence  the  chosen  S,-  and  Sjm  step  2e  of 
the  algorithm  satisfy  ISj^jSy  I S(3/4)- IP  I .[] 

Step  1  of  the  algorithm  takes  G(n-fm)  time  and  renders  G  sparse  (Pin)  arcs).  As  mentioned 
above,  steps  2a  and  2b  can  be  implemented  to  run  in  G(/j)  time  usirtg  the  algorithm  in  [Ta2],  in 
conjunction  with  the  algorithms  in  [HT]  and  [GT].  Each  of  steps  2c  through  f  takes  0(n)  time. 
Hence  each  execution  of  the  repeat  loop  takes  linear  time.  Since  by  Lemma  1 1  the  repeat  loop  is 
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executed  OOogn)  times,  the  enure  transitive  reduction  algorithm  runs  in  0(m+nlogrt)  time. 

The  algorithm  of  section  2  can  also  be  implemented  to  run  in  Oim+n\ogn)  time.  This  is 
because  the  minimum-weight  branching  algorithm  of  Edmonds  [Ed2]  can  be  implemented  to  mn 
in  linear  time  for  0-1  edge  weights  by  using  the  algorithm  in  [GGST],  with  the  heaps  replaced  by 
two  buckets.  As  before,  the  redundant  arcs  can  be  found  in  linear  time  and  hence  each  execution 
of  the  repeat  loop  takes  linear  time,  leading  to  an  O  (m+nlogn)  time  sequential  algorithm  for  tran¬ 
sitive  reduction. 

We  have  obtained  sequential  and  parallel  algorithms  with  similar  complexities  for  analo¬ 
gous  problems  on  undirected  graphs,  i.e.,  for  finding  a  minimal  bridge-connected  spanning  sub¬ 
graph  and  a  minimal  biconnected  spanning  subgraph  in  an  undirected  graph,  if  such  subgraphs 
exist  These  results  will  appear  in  a  companion  paper. 

We  conclude  by  noting  that  it  is  conceivable  that  one  (or  both)  of  our  sequential  algorithms 
mns  in  linear  time,  since  it  is  possible  that  the  repeat  loop  needs  to  be  executed  only  a  constant 
number  of  times.  We  leave  this  question  for  further  investigation.  For  the  same  reason  it  is  possi¬ 
ble  that  our  parallel  algorithms  tun  faster  than  the  stated  time  bounds  by  an  O  Gogn)  factor. 
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